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ABSTRACT 

Understanding the stability of the magnetic field in radiation zones is of cru- 
cial importance for various processes in stellar interior like mixing, circulation 
and angular momentum transport. The stability properties of a star containing 
a prominent toroidal field in a radiation zone is investigated by means of a linear 
stability analysis in the Boussinesq approximation taking into account the effect 
of thermal conductivity. The growth rate of the instability is explicitly calculated 
and the effects of stable stratification and heat transport are discussed in detail. 
It is argued that the stabilizing influence of gravity can never entirely suppress 
the instability caused by electric currents in radiation zones although the stable 
stratification can significantly decrease the growth rate of instability 

Subject headings: MHD - instabilities - stars: interiors - stars: magnetic fields - 
Sun: magnetic fields 

1. Introduction 



Magnetic fields localized in stellar radiation zones can play an important role in many 
essential phenomena like mixing, angular momentum transport and formation of the solar 
tachocline, for instance. 
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It is rather uncertain which field can be present in the radiation zone. Observations 
provide upper limits on the allowed strength of the magnetic field in the radiation zone of 
the Sun. For instance, helioseismological measurements suggest an upper limit B < 4 x 10 7 
G as an order of magnitude estimate (see, e.g., Friedland & Gruzinov 2004). The recently 
measured oblateness of the Sun (Codier & Rozelot 2000) implies that the strength of the 
magnetic field is lower than 7 x 10 6 G. The observations of the splitting of solar oscillation 
frequencies provide the upper limit of B < 3 x 10 5 G for a toroidal field near the base of 
the convective zone (Antia et al. 2000). For other stars estimates are less certain and only 
theoretical upper limits can be derived. 

The possible origin of the field (if it exists) is also unclear. Likely, a hydromagnetic 
dynamo cannot operate in the radiation zone, where no strong flows are available to sustain 
a vigorous dynamo action. Perhaps, relic magnetic fields acquired by the star at the early 
stage of evolution can persist there. These type of fields could have formed, for instance 
because of differential rotation, which could have stretched a weak diffuse primordial seed 
field (see, e.g., Dicke 1979) into a dominant toroidal field. Due to a high conductivity, the 
ohmic decay is very slow and the decay time can exceed several billion years. Therefore, once 
formed, the large-scale relic field would survive in the radiation zone during the life-time of 

db SlclX. 

The magnetic field, however, can evolve in a radiation zone not only due to ohmic 
dissipation but also because of the development of various instabilities. For instance, in 
differentially rotating radiation zones, the magnetorotational instability can occur. Most 
likely, however, the magnetized radiative zones rotate rigidly and the stability of magnetic 
configurations is determined by various current-driven instabilities. Such instabilities are 
well studied in cylindrical geometry in the context of laboratory fusion research (see, e.g., 
Freidberg 1973, Goedboed 1971, Goedbloed & Hagebeuk 1972). For example, the stability 
properties of a pure toroidal field B v are determined by the parameter a = d In B v (s)/d Ins 
where s is the cylindrical radius. The field is unstable to axisymmetric perturbations if a > 1 
and to non-axisymmetric perturbations if a > —1/2 (Tayler 1973a,b, 1980). Note, however, 
that the addition of a even relatively weak poloidal field alters substantially the stability 
properties of the magnetic configuration. For example, if the poloidal field is uniform and 
relatively weak, the instability condition reads a > — 1 at variance with the condition of 
instability for a purely toroidal field (see, e.g., Knobloch 1992; Dubrulle & Knobloch 1993) 
which predicts that an unstable toroidal field configuration has a > — 1/2. In astrophysical 
conditions, the instability caused by electric currents might have various characteristic prop- 
erties (see Bonanno & Urpin 2008a,b). In the presence of both azimuthal and axial fields, 
non-axisymmetric disturbances with large azimuthal wavenumbers m turn out to be most 
rapidly growing. Unstable disturbances exhibit a resonant character, i.e. the wave vector 
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k = (m/s)e<p + k z e z approximately satisfies the condition of magnetic resonance, B ■ k = 

— * 

where k z is the wavevector in the axial direction and B is the magnetic field. The length 
scale of this instability depends on the ratio of poloidal and azimuthal field components and 
it can be very short, while the width of the resonance turns out to be very narrow. For this 
reason, its excitation in direct numerical simulations can be problematic. 

Stability of the spherical magnetic configurations is studied in much less detail and even 
the overall stability properties of radiation zones are rather unclear. Braithwaite & Nordlund 
(2006) have studied the stability of a random initial field in the stellar radiative zone by means 
of direct numerical simulations and found that the stable magnetic configurations generally 
have the form of tori with comparable poloidal and toroidal field strengths. Numerical 
modeling by Braithwaite (2006) confirmed that the toroidal field with oc s or oc s 2 is 
unstable to the m — 1 mode as it was predicted by Tayler (1973). However, even a purely 
toroidal field can be stable in the region where it decreases rapidly with s. Note that a purely 
toroidal field cannot be stable throughout the whole star because the stability condition for 
axisymmetric modes (a < 1) is incompatible with the condition that the electric current 
in the z-direction has no singularity at s — > which implies a > 1. The stability of the 
toroidal field in rotating stars has been considered by Kitchatinov (2008) and Kitchatinov 
& Riidiger (2008) who argued that the magnetic instability is essentially three-dimensional 
and determined the threshold field strength at which the instability sets. Estimating this 
threshold in the solar radiation zone, the authors impose the upper limit on the magnetic 
field « 600 G. 

In this paper, we consider the stability of the toroidal field in radiation zones by taking 
into account stratification and thermal conductivity. It turns out that magnetic configu- 
rations can be stable or unstable depending on the radial profile of the toroidal field. We 
argue that stable stratification can substantially decrease the growth rate of the instability 
but cannot suppress it entirely. 



2. Basic equations 

Consider the stability of an axisymmetric toroidal magnetic field in the radiation zone 
using a high conductivity limit. We work in spherical coordinates (r, 9, ip) with the unit 
vectors (e r , e#, e^). We assume that the toroidal field depends on r and 9, B v = B v (r,9). 
In the incompressible limit, the MHD equations read 

^ + (v- V)v = + g + -!-(V xB)xB, (1) 
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dB 



dt 



-Vx(nB)=0, 



(2) 

V-v = 0, V-B = 0, (3) 
where g is gravity. In the basic state, the gas is assumed to be in hydrostatic equilibrium, 



then 



Vp 1 - 

— =g + (V x B) x B. 

p Anp 



(4) 



We assume that the magnetic energy is subthermal and therefore g is approximately radial. 
The equation of thermal balance reads in the Boussinesq approximation 



^ + t?-(VT- V ad T) = V • (kVT), 
where k is the thermal diffusivity and V a dT is the adiabatic temperature gradient. 



(5) 



We consider a linear stability. Small perturbations will be indicated by subscript 1, 
while unperturbed quantities will have no subscript. Linearizing Eqs.(l)-(3) and (5), we 
take into account that small perturbations of the density and temperature in the Boussinesq 
approximation are related by pi/p = — /3(Ti/T) where /3 is the thermal expansion coefficient. 
For small perturbations, we use a local approximation in the ^-direction and assume that their 
dependence on 9 is proportional to exp(— U9), where I ^> 1 is the longitudinal wavenumber. 
Since the basic state is stationary and axisymmetric, the dependence of perturbations on t 
and ip can be taken in the exponential form as well. Then, perturbations are proportional to 
exp (at — U9 — imip) where m is the azimuthal wavenumber. The corresponding wavevectors 
are kg = l/r and k v = m/r sin 9, respectively. The dependence on r should be determined 
from Eqs.(l)-(3), (5). For the sake of simplicity, we assume that unperturbed p and T are 
approximately homogeneous in the radiation zone. This assumption does not change the 
main conclusions qualitatively but substantially simplifies calculations. We also neglect the 
effect of rotation. As we know from the work of Pitts and Tayler (1986), rotation itself cannot 
remove instabilities of the interchange type, although it can affect the shape the unstable 
modes. For this reasons we argue that the conclusions of our investigation cannot be altered 
by the inclusion of rotation. 

Eliminating all variables in favor of v\ r and T±, we obtain the following set of two coupled 
equations 

" 2 



(a 2 + c4K + {^a 2 + ^u>£j v' lr + 



y-kl(v 2 + u; 2 A ) 



(6) 



2 ^ { \ }z \ 



2 k 2 a 2 



A/ J 



+ r UA \Hkl rk^a 2 ^ ,; 2 



vir = -kiPga-^, 



T 



3 

T 



U BV 



-Vir, 



(7) 
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where the prime denotes a derivative with respect to r and 

47TP' UBV ~ T 



^ = ^, ujl v = - 9 4(V ad T-VT) r , (8) 



2 2 2 1 d 
k ± — k 9 + ky, — = —(rB v ). 

Some general stability properties can be derived directly from Eqs.(6)-(7). Consider 
perturbations with a very short radial wavelength for which one can use a local approximation 
in the radial direction, such as v lr oc exp(— ik r r), where k r is the radial wavevector. If 
k r ^> max(A; , k v ), then Eqs.(6)-(7) can be reduced with the accuracy in terms of the lowest 
order in (/c r r) _1 to the following set 

-(* + ^ = ^v lr , A; r V + u%)v lr = ki(3ga^, (9) 
where k 2 = k 2 + k\. The corresponding dispersion relation reads 

„2_2 , / . .2 , h ±, .2 \ , t2, ,2 



cr + kA;V 2 + ^ + -j^^Bvj <* + Kk 2 u 2 A = 0. (10) 

The conditions that at least one of the roots has a positive real part (unstable mode) is 
determined by the Routh criterion (see, e.g., Aleksandrov, Kolmogorov, & Laurentiev 1985). 
For a particular case of a cubic equation, these conditions can be easily obtained, for exam- 
ple, from expressions derived by Urpin & Riidiger (2005). Since the quantities k and u\ are 
positively defined, the only non-trivial condition of instability is oo 2 BV < that is not sat- 
isfied in the radiation zone by definition. Therefore, modes a with short radial wavelength 
are always stable to the current-driven instability contrary to the conclusion obtained by 
Kichatinov (2008) and Kichatinov & Riidiger (2008). 



3. Numerical results 

We assume that the radiation zone is located at Ri > r > R. The toroidal field can be 
represented as 

B ip = B il)(x)sm8, (11) 

where x = r/R and Bo is the characteristic field strength; ip ~ 1 is a function of the radius 
alone. The dependence of ip on x is uncertain in in the radiation zone and, in this work, we 
consider only the case where ip increases with x. Other possibilities will be considered in a 
forthcoming paper. We parametrize ip(x) with the following dependence 

i>(x) = (d x exp + £> 2 ) (12) 
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where d characterizes the width of a field distribution and D\ and D 2 are chosen so that 
■ip(x = 1) = 1 and i]){x = Xj) = as shown in Fig.l; x^ = Ri/R. We choose Xj to be equal to 
0.01 from computational reasons. We have verified that our results are basically insensitive to 
the precise value of Xj as long as it is close to the center. The situation where the field reaches 
its maximum at the outer boundary can model, for example, the radiative interior of a star 
with a convective envelope. In this case, the bottom of the convection zone is the location of 
the toroidal field generated by a dynamo action which can penetrate into the radiation zone. 
This situation can also mimic the toroidal field in the liquid core of neutron stars. Likely, the 
magnetic field of these objects is generated by turbulent dynamo during the very early phase 
of the evolution when the neutron star is subject to hydrodynamical instabilities (see, e.g., 
Bonanno et al. 2003). Dynamo induced by turbulent motions generates the magnetic field 
of complex topology including small scale fields (Urpin & Gil 2004). Large-scale dynamo is 
most efficient in the surface layers where the density gradient is maximal. Therefore, the 
generated field increases outward and reaches its maximum in the outer layers (Bonanno et 
al. 2005, 2006). This magnetic field can be subject to current-driven instabilities after the 
end of the unstable phase. 

Introducing the dimensionless quantities 



4irpR% ujao u A o ^ao 



where lot = k/R%, we can transform Eqs.(6)-(7) into a dimensionless form. These equations 
with the corresponding boundary conditions describe the stability problem as a non-linear 
eigenvalue problem. Fortunately, the main qualitative features of this problem are not sen- 
sitive to the choice of boundary conditions. That is why we choose the simplest conditions 
and assume that v\ r = T\ = at r = Ri and r = R. Note that the parameter 5 is large in 
radiation zones but, most likely, e is relatively small if the magnetic field is not very weak. 
In calculations, we suppose 5 and e to be constant through the radiation zone. 

To illustrate the stabilizing effect of stratification on the instability, we plot in Fig. 2 
the growth rate as a function of the stratification parameter 5 for several eigenmodes with a 
different number of nodes in the radial direction. In these calculations, we neglect the thermal 
conductivity, therefore, the stabilizing effect of gravity is most pronounced. The growth rate 
is always maximal for the fundamental eigenmode (n = 0) but rapidly decreases with an 
increasing number of nodes n. Note that such a rapid decrease of Y with the number of an 
eigenmode is typical also for the case of neutral stratification 5 = (or no gravity g — 0). 
The conclusion that T decreases rapidly with n confirms our result obtained from Eqs.(6)- 
(7) in a short wavelength approximation (see Section 2) that modes with a short radial 
lengthscale (large n) should be stable. This conclusion is at variance with the statement by 



-7- 



Kitchatinov & Riidiger (2008) that the most rapidly growing modes correspond ton~ 10 3 . 
According to our calculations, the fundamental mode turns out to be most rapidly growing. 
However, even the instability of this mode is entirely suppressed if 8 > 9. The effect of a 
thermal conductivity can change the properties of a current- driven instability qualitatively. 
In Fig.3, we plot the dependence of the growth rate on the parameter stratification 5 2 for 
three different values of the thermal conductivity, corresponding to e — 10~ 2 , 2 x 10~ 2 , and 
4 x 10~ 2 . The behaviour of all curves is qualitatively similar: the growth rate is « 1 at small 
8 decreases as 

T oc (T 2 (14) 

for 5 2 > 100 or, in dimensional form, 

a oc u A o(uj ao /ubv) 2 - (15) 

This dependence can be obtained also directly from the basic equations. The stabilizing 
effect of gravity in a linear theory is determined by the first term on the right hand side 
of Eq.(l). Linearization of this term yields (Vp/p)i « —<JPi/p ~ gTi/T since pi in 
the Boussinesq approximation. Perturbations of the temperature are determined by thermal 
balance Eq.(7). Near the threshold of the instability when a is small, we have T\jT « 
— (u 2 3V /(3gK,k 2 _)vi r . Therefore, the stabilizing contribution of gravity in Eq.(l) is of the 
order of (oj^v / K ^±) v ^r ~ 8 2 (u\ / K,k"j_)v lr . On the other hand, the destabilizing effect of 
electric currents in Eq.(l) is given by the Lorentz force. The order of magnitude estimate 
of the Lorentz force yields ~ B^Bi^/Anpr. Since the toroidal field is inhomogeneous in the 
basic state, perturbations of the toroidal field are produced basically by perturbations of the 
radial velocity and B lip ~ {B^jaH^vxr where H is the radial lengthscale of the magnetic field. 
Using this expression, we obtain that the Lorentz force is of the order of ~ (uj 2 m / aa)v\ r where 
a = H/R. Equating the stabilizing contribution of gravity to the destabilizing contribution 
of the electric current, we obtain 

*~U A o ^ (k 2 ± R 2 2 ). (16) 

If e is fixed, this expression yields the dependence Y oc 5~ 2 that follows from numerical 
calculations. 

The relation (14) implies that even a very strong stable stratification cannot entirely 
suppress the instability. The growth rate turns out to be non-vanishing even for large 5 at 
variance with the case where thermal conductivity is neglected (see Fig. 2). 
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4. Conclusions 

We have considered the effects of stratification and thermal conductivity on the stability 
of a toroidal field in stellar radiation zones. The magnetic configuration with a predominantly 
toroidal field can be formed, for example, due to differential rotation during the early stage 
of stellar evolution if the star has got even a weak seed field. 

It turns out that stable stratification can suppress the current induced instability of the 
toroidal field if perturbations are not influenced by the thermal conductivity (very small e). 
The instability does not arise if the Brunt- Vaisala frequency is greater than ~ 9uao- Since 
ubv is typically high in radiative zones (~ 10~ 3 — 10~ 4 s _1 ) the instability sets in only if 
the field is very strong (> 10 6 — 10 7 G). Higher eigenmodes are suppressed stronger than the 
fundamental one and perturbations with short radial wavelength are always stable. 

The thermal conductivity drastically changes the character of the instability. This 
concerns particularly the behavior near the threshold of instability where the growth rate 
is small. It turns out that the growth rate is non-vanishing for any stratification. Even 
very strong gravity cannot suppress the instability entirely but it only decreases the growth 
rate. This sort of behavior is typical for perturbations with any wavevector k± and can 
be easily understood. A destabilizing effect of electric currents in the momentum equation 
(1) originates from the last term on the r.h.s. and is proportional to a perturbation of the 
magnetic field, B\. A perturbation B\ is related to v\ by Bi a tii/tx as it follows from the 
linearized induction equation. A stabilizing influence of gravity in the linearized momentum 
equation is proportional to a perturbation of the temperature, 7\. If the growth rate is 
small, then we can obtain from thermal balance equation (5) that 7\ oc vi/u T - Comparing 
the destabilizing and stabilizing effects, we see that stable stratification can never suppress 
the instability at variance with the k = case. 

A decrease of a caused by stratification is inversely proportional to the Brunt- Vaisala 
frequency. If gravity is strong but the magnetic field is weak, the instability develops very 
slowly. Generally, for a sufficiently weak magnetic field, the growth rate can be comparable 
to the inverse life-time of a star. It should be noted also that, most likely, the field does not 
decay to zero because of this instability. When the field becomes weaker, the growth rate 
of the instability decreases and the field cannot decay to values smaller than those resulting 
from a growth rate of the order of the inverse life-time of a star. Therefore, a weak field can 
generally be only slowly changed during the life of the star, although its radial profile can 
be unstable. 

VU thanks also INAF-Ossevatorio Astrofisico di Catania for hospitality and financial sup- 
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Fig. 1. — The dependence of the toroidal field on the spherical radius for d = 0.2 (solid line), 
0.1 (dashed), and 0.05 (dot-dashed). 
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Fig. 2. — The dimensionless growth rate as a function of 5 for several eigenmodes with 
various number of nodes n in the radial direction, for m — 1, I — 20, d — 0.1 and e — 0. 
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Fig. 3. — The dimensionless growth rate as a function of 5 2 for the fundamental eigenmode 
with m = 1, I = 20, d = 0.1, and three values of e: 10~ 2 (solid line), 2 x 10~ 2 (dashed), and 
4 x 10~ 2 (dot-dashed). 



